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Abstract 

In this work we provide a novel approach for computing the coefficients 
of the characteristic polynomial of a square matrix. We demonstrate that 
each coefficient can be efficiently represented by a set of circle graphs. Thus, 
one can employ a diagrammatic approach to determine the coefficients of the 
characteristic polynomial. 

1 Introduction 

A variety of different branches in Mathematical Physics boil down to calculating the eigenvalues 
of an n X n matrix A over some field K usually taken to be the real line or the complex plane. 
When we are interested in all the latent roots of the characteristic polynomial of a matrix and 
not in its latent vectors we usually expand out the expression [1] 

n 

/(x) = det{A- xl) = Y[{x - Xi) 

i=l 

= + (-l)^ei(A)x"-i + (-l)2e2(A)x"-2 + {-lfe3{\)x''-'' + ■■■ + (-l)"e„(A) 

(1) 

where efc(A) = J2l<^l<■■■<^^.<n ' ' ' ^ik elementary symmetric polynomials in k eigen- 

values Ai^, s = l,---,k. There are different methods to determine explicitly efc(A) in terms 
of a sum of trace products nt=i itr{A^))"^'' , over all possible partitions • ■ • , Z™'') of the 
positive integer k. Undoubtably, all methods for high order matrices become very tedious. 
One method for an n x n matrix is based upon the multinomial formula 

(Ai + --- + A„r= E f . . (2) 

where the summation is taken over all sequences of nonnegative integer indices ki, ■ ■ ■ ,kn such 
that J27=i ki = m and the cofficients are given by 

( ]= (3) 

\k^,---,kn J ki\---kn\' ^ ^ 

^We adopt the following notation: n = 2'"^, • • • , r™"^, • ■ •) where nii counts the number of parts of n 

which are equal to i. It is called the multiplicity of i in n. 
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For n = 2 relation ([2]) gives 

E^n) = E An A.. (4) 

from which we find 

e,{X) = -U{tr{A') - (trAr)) . (5) 



Similarly for 63 (A) we have 



\ 



\ n n 

EAn -2EAf,-3 E A^^A, 

il=l / «i=l 1^^42 = 1 



-i (2tr{A^) + 3tr(A2)trA - (trA)^) . (6) 
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Note that the three terms in (jH]) correspond to the partitions 3^ = (1^,2^) = 1^. In the same 
vein we can compute any efc(A). 

A second method is to use Newton — Girard's formulas which give the connection between 
the coefficients Pr-(A) = (— l)~^er(A) and the power sums 



'S'r - E Ail, 
n=l 

Sr + Sr-lPl+ ■ ■ ■ +SiPr-l + = 0, T < U. (7) 

From it follows that 

Pi = -Si 

P2 = -^iSipi + S2) 

Pr = --(S'lPr-lH \- Sr^iPi + Sr). (8) 

r 

Relations ([H]) determine the coefficients pr and the process is called Le Verrier's method. 



2 The New Method 

The method we propose for the evaluation of the coefficients of the charactristc polynomial is 
based upon the knowledge of the number of partitions q{k) for the index k which specifies the 
elementary symmetric polynomial efc(A). The generating function of q{k) is given by Euler's 
formula [2] 

-1 oo 

Ftotix) = ... =EgW^^- (9) 



Other useful generating functions are 

1 , _ 

k=l 



I I fc 

Feven{x) = j-joo ^^-^ _ ^2k^ ' ^oddix) = j-joo _ ^2k+l^ ' ^uneq.ix) = + X ) (10) 
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where F^ven-, Fodd and Fumq. enumerate partitions of k into even, odd and unequal parts respec- 
tively. 

Having the partitions at hand we proceed by imposing certain rules for the construction of 
diagrams: 

(1) Each matrix (Aij) is represented by a line segment with indices i,j attached to the 
endpoints. The trace of A is graphically formed by gluing the two endpoints and thus 
resulting in a circle graph. 

(Aj) i j 

tT{A) O 

(2) For every positive integer k we associate a single circle graph with k points (we call it a 
k-circle from now on) and contributing a factor —tr{A^). 

(3) A partition of A; = (1™% 2"^^, • ■ • , r"*'', • ■ •) is represented by a set of circle graphs which 
can be constructed by a cutting and sewing procedure from the k-circle 

O -(orHo)-v..(o)-v.. 

fc— circle 

where the powers on the right handside stand for the multiplicity of each graph. 

(4) If r™'' is a single partition of k then the coefficient is given by 

The interpretation of each factor in ( fTTl) is as follows: 

(a) The number of circular permutations of r distinct points on the circle is (r — 1)!. 

The number of permutations of r-circles nir times is (m^)!. 

(7) The ways of extracting r points each time from a fc-circle, times, without replacement 
and disregarding order is given by the product. 

Note that if we sum up the absolute values of factors (fTT]) for all possible partitions of k then 
we recover k\. 

A more involved case study will be the partition k = (r'"'",s™^). The combinatorial factor 
now reads 

i-^r--^i-MA, ^'-ZZ^r' E n ( T - ) ■ (-) 

As an application consider the case of k = 6. The total number of partitions is q{6) = 11 
which splits into QeveniG) = 3 and QoddiG) = QuneqX^) = 4 parts. Adopting the convention 
T;'^ = (tr(A'))^ a list of the contributions of all partitions of 6 is given in the following table 
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Table 1: The summands of cq applying the diagrammatic approach. 



Partitions 


Coefficients 


(6) 


-120 Tg 


(5,1) 


144 T5T1 


(4,2) 


90 T4T2 


(32) 


40 


(4,12) 


-90 T^Tf 


(3,2,1) 


-120 TgTaTi 


(23) 


-15 T| 


(3,1^) 


40 TsTf 


(2M2) 


45 T^T^ 


(2,1^) 


-15 T2Tf 


(1^) 





The connection of these coefficients with the elementary symmetric polynomial eg (A) is 

p,{X) = (-i)%(A)= 

l<il<---<i6<n 

= (120 Tq - 144 T5T1 - 90 TiT2 - 40 r| + 90 TJ'l + 120 TsTsTi + 15 
6! \ 

-40 T^Tf - 45 r|T2 + 1 5 TsT^^ - Tf ^ ( 13) 

3 Conclusions 

In this letter we present explicitly a diagrammatic way to calculate the coefficients of the 
characteristic polynomial of a square matrix. All information is encoded in combinatorial form 
into the sets of circle graphs constructed for all partitions of the index associated with the 
corresponding elemetary symmetric polynomial. 
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